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CRRESPRO  Documentation 


1.  OVERVIEW 

CRRESPRO  predicts  proton  omnidirectional  fluence  per  year  and  integral  omnidirectional  fluence  per 
year  at  selected  energies  in  the  range  1  •  100  MeV  for  an  orbit  specified  by  the  user.  It  closely  parallels 
its  counterpart,  CRRESRAD,  which  predicts  dose  behind  four  thicknesses  of  hemispherical  aluminum 
shielding.  Likewise,  this  documentation  closely  follows  the  CRRESRAD  documentation.  The 
CRRESPRO  software  uses  flux  models  created  from  data  collected  by  the  proton  telescope  (PROTEL) 
on  board  the  Combined  Release  and  Radiation  Effects  Satellite  (CRRES)  flown  from  25  July  1990  to  12 
October  1991  during  solar  maximum.  CRRES  was  in  a  geosynchronous  transfer  orbit  with  an  inclination 
of  18°,  a  perigee  of  350  km,  and  an  apogee  of  33000  km.  It  traversed  the  radiation  belts  twice  per  orbit 
with  a  period  of  9  hours  52  minutes.  In  March  1991,  a  magnetic  storm  caused  a  reconfiguration  of  the 
inner  magnetosphere,  resulting  in,  among  other  features,  double  proton  h<  its  forming  over  a  certain 
energy  range.  Because  of  this  change,  two  CRRES  models  were  created.  The  quiet  model  uses  data 
from  July  1990  to  March  1991,  and  the  active  model  uses  data  from  March  1991  to  October  1991.  Note 
that  in  this  documentation  and  the  CRRESPRO  software,  "quiet"  refers  to  the  period  from  July  1990  to 
March  1991  (single  proton  belt)  and  "active"  refers  to  March  1991  to  Octobet  1991  (double  proton  belt). 
Quiet  and  active  as  used  here  for  the  inner  radiation  belt  have  no  correspondence  to  quiet  and  active  as 
determined  by  Kp.  In  fact,  the  average  Kp  for  the  two  CRRES  periods  was  the  same,  namely  2.2. 

The  CRRESRAD  software,  which  was  released  by  this  office  (PL/GPSP)  in  1992,  was  developed 
when  dose  measurements  on  CRRES  and  the  Defense  Meteorological  Satellite  Program  (DMSP)  showed 
that  the  current  NASA  models  used  to  predict  dose  for  satellite  design  and  orbit  determination  are  not 
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always  adequate  (Oussenhoven  et  al. ,  1987,  1991a,  1991b).  A  major  geomagnetic  storm  occurred  in  late 
March  1991  following  an  intense  solar  proton  event  during  which  strong  shock  fronts  hit  the 
magnetosphere.  The  result  of  the  shock  passages  was  a  change  in  the  magnetosphere  to  a  more  active 
configuration  (Mullen  et  al.,  1991  and  Blake  et  al.,  1992).  Before  the  March  storm,  dose  predicted  by 
the  NASA  models  for  protons  below  100  MeV  agreed  very  well  with  measured  dose  in  the  inner  zone, 
but  dose  predicted  by  the  NASA  models  for  electrons  was  larger  than  the  measured  dose  by  as  much  as 
an  order  of  magnitude.  After  the  storm,  a  second  energetic  proton  belt  and  a  third  energetic  electron  belt 
were  formed  in  the  slot  region  where  there  had  previously  been  little  radiation.  During  this  period,  mqjor 
differences  exist  between  the  CRRES  and  NASA  models  (Brautigam,  1992;  Oussenhoven,  et  al.,  1992). 

A  well  designed,  well  calibrated  proton  detector  can  be  used  to  create  accurate,  pitch  angle 
dependent,  empirical  proton  flux  models  on  the  magnetic  equator  as  a  function  of  L  shell  (Bee  below)  for 
a  wide  energy  range.  These  average  fluxes  can  then  be  propagated  down  the  relatively  stable  magnetic 
field  lines  of  the  inner  magnetosphere  to  extend  the  models  outside  the  range  of  the  satellite.  A  utility 
such  as  CRRESPRO  will  give  the  total  proton  fluence  predicted  by  the  model  at  given  energies  for  an 
arbitrary  satellite  orbit.  Then,  with  a  transport  code,  the  results  of  the  model  can  be  used  to  determine 
dose  behind  any  type  of  shield.  Furthermore,  different  models,  such  as  the  CRRES  proton  models  and 
the  NASA  proton  models,  can  be  directly  compared  for  differences  in  intensity,  pitch  angle  dependence, 
and  spectral  shape. 

Such  a  comparison  was  made  by  Gussenhoven,  et  al.  Among  the  findings  of  the  comparison  are 
two  regions  in  the  inner  magnetosphere  that  are  populated  by  high  energy  protons.  The  inner,  stable 
region  exists  below  L  -  1.8.  Beyond  L  ■  1 .8,  there  exists  a  population  that  was  previously  thought  to 
be  a  continuous  extrapolation  of  the  inner  population,  as  is  shown  by  the  NASA  AP8MAX  model. 
However,  the  PROTEL  measurements  show  that  beyond  L  -  1 .8  a  new  population  can  be  brought  in 
during  periods  of  extreme  solar  wind  conditions  and  apparently  added  to  the  existing  population.  This 
presents  a  radiation  hazard  that  can  continue  for  months  and  is  not  dealt  with  by  the  NASA  models.  It 
is  possible  that  AP8MAX  is  an  average  of  the  two  conditions  measured  by  CRRES.  (Gussenhoven  et  al. , 
1993) 

CRRESPRO  requires  an  80386  or  80486  DOS  based  computet  with  a  math  coprocessor,  a  color 
VGA  card,  and  a  hard  drive.  A  printer  is  not  needed,  but  it  is  useful.  The  active  directory  must  be  the 
one  containing  CRRESPRO.EXE  and  its  suppor  files,  CRRESPRO  is  run  by  typing  "CRRESPRO"  on 
the  command  line.  The  user  will  be  queried  for  the  name  of  the  orbit  to  be  evaluated.  This  name  will 
be  used  for  all  of  the  data  files  subsequently  created  for  this  orbit.  The  data  files  can  be  written  to 
another  directory  by  adding  the  path  to  the  name,  For  example,  to  run  an  orbit  called  "Sample"  and  store 
the  data  In  "C:\Orbits\",  the  name  would  be  specified  as  "C:\Orbits\Satnple".  No  extension  can  be  given 
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in  the  name  as  default  extensions  are  used.  An  ASCII  README  file  is  included  with  the  CRRESPRO 
software.  To  read  this  file  enter  "TYPE  README  |  MORE”  on  the  command  line  and  press 
<  enter  >.  README  documents  recent  updates  to  the  software  and  has  a  directory  list  of  the  files 
included  with  CRRESPRO.  As  a  precaution  against  viral  contamination,  make  sure  the  ".EXE"  files  are 
the  same  size  and  date  as  given  in  the  README  file.  A  sample  session  for  determining  the  fluence  from 
a  CRRES-like  orbit  is  shown  in  Appendix  A. 


2.  FLUX  MODEL 

The  motion  of  the  high  energy  ions  and  electrons  that  cause  dose  deposition  in  satellites  in  Earth’s 
magnetosphere  is  determined  almost  entirely  by  Earth’s  magnetic  field.  These  particles  are  trapped  by 
the  magnetic  field,  and  are  constrained  to  bounce  back  and  forth  between  the  north  and  south  poles  along 
magnetic  field  lines.  The  paths  these  particles  take  are  referred  to  as  "shells."  Because  these  particles 
are  constrained  to  move  along  magnetic  field  lines,  fluence  resulting  from  these  particles  is  best  ordered 
by  specifying  position  with  magnetic  field  parameters  instead  of  using  spatial  coordinates.  The  magnetic 
field  parameters  used  for  the  flux  model  are  the  Mcllwain  L  parameter  (Mcllwain,  1961),  which  is 
referred  to  hereafter  as  L  shell,  and  B/Bq.  In  a  dipole  magnetic  field  the  L  shell  of  a  given  point  in  space 
is  the  equatorial  radial  distance  in  Earth  radii  (Rg)  of  the  field  line  passing  through  the  given  point. 
Earth’s  magnetic  field  has  some  non-dipole  components,  and  L  shell  is  calculated  by  determining  which 
dipole  field  line  the  actual  field  line  would  correspond  to  if  Earth's  magnetic  field  were  relaxed  to  a 
dipole  field.  B/B0  is  the  ratio  of  the  magnetic  field  magnitude  at  a  point  in  space  divided  by  the  minimum 
magnitude  of  the  magnetic  field  on  the  same  field  line.  The  field  is  strongest  near  the  poles  and  weakest 
near  the  equator  so  that  B/B0  »  1.0  is  near  the  equator  and  maximum  B/B0  is  closer  to  the  poles.  For 
a  dipole  field,  B/B„  is  a  function  only  of  magnetic  latitude.  The  solid  lines  in  Figure  1  show  L  shells  2, 
3, 4.  3,  and  6  of  a  dipole  magnetic  field  centered  in  Earth.  The  dotted  lines  show  constant  B/B„  contours 
and  their  magnitude  spaced  every  10°  in  an  Earth-centered  dipole  field  between  -40°  and  +40°  magnetic 
latitude.  Chapter  5  of  the  Handbook  of  Geophysics  and  the  Space  Environment  (Spjeldvik  and  Rothwell, 
1983)  contains  a  more  complete  description  of  Earth's  magnetosphere  and  radiation  belts. 

The  Olson-Pfitzer  static  external  magnetic  field  model  and  the  IGRF85  internal  magnetic  field  model 
were  uced  to  calculate  the  L  shell  and  B/Bu  values  when  binning  the  CRRES  flux  data.  The  IGRF83 
model  specifies  the  magnetic  field  that  comes  from  currents  inside  Earth  (IAGA  Division  I,  Working 
Group  1,  1986).  This  field  is  nearly  dipolar  with  an  1 1.3s  tilt  and  436  km  offset  with  respect  to  Earth’s 
spin  axis  and  is  slowly  changing  vith  time  (Knrcht  and  Shuman,  1985).  The  IGRF  model  accounts  for 
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Figure  1.  L  Shells  2,  3,  4,  5,  and  6  (Solid  Lines)  and  Lines  of  Constant  B/Bt  -  1.0,  1.1,  1.7,  3.1,  and  7.4 
(Dashed  Lines)  in  a  Dipole  Field  Centered  in  Barth. 

the  alow  time  variation  of  the  field  as  well  as  the  tilt  and  non-dipole  contributions  to  the  field.  The 
Olson-Pfitzer  static  model  (Olson  and  Pfltzer,  1974)  predicts  changes  in  the  magnetic  field  caused  by 
currents  in  the  magnetopause,  the  ring  current,  and  the  magnetosphere's  tail.  The  model  also  recounts 
for  the  tilt  of  Earth's  dipole  with  respect  to  the  Earth-Sun  line.  The  external  currents  that  contribute  to 
the  magnetic  field  change  with  variations  in  the  solar  wind  and  magnetospheric  activity.  The  static  model 
assumes  these  currents  do  not  change  and  uses  values  from  average  quiet  conditions.  We  should  note 
that  the  Olson-Pfitzer  external  field  makes  no  contribution  inside  an  L  shell  of  2.4. 

The  flux  models  used  by  CRRESPRO  are  based  on  in  situ  flux  measurements  made  by  CRRES. 
Their  construction  is  described  in  Gussenhoven  et  al  (1993),  and  is  briefly  reviewed  here.  The 
Instrument  used  to  measure  flux  on  CRRES  was  the  Proton  Telescope  (PROTEL).  PROTEL  had  two 
detector  heads,  which  together  measured  protons  from  1  to  100  MeV  in  24  energy  steps,  giving  a 
complete  spectrum  every  1 .024  seconds.  The  angular  resolution  of  the  detector  low  (high)  energy  head 
was  ±  10°  x  ±10°  (±  12°  x  ±  17°).  A  ftill  description  of  PROTEL  is  available  in  Violet  et  al.  (1993) 
and  Lynch  et  al.  (1989).  The  PROTEL  detectors  comprised  detector  stacks  and  a  logic  system  that 
required  single  or  double  coincidence  to  verify  that  the  proper  energy  particle  is  counted.  In  addition, 
both  active  and  passive  shielding  were  used  around  much  of  the  detector  stack.  The  detectors  were 
extensively  calibrated  prior  to  launch.  During  calibration,  we  found  that  energetic  protons  ( >60  MeV) 
incident  over  a  large  angular  range  with  respect  to  the  detector  axis  could  degrade  sufficiently  In  the 
shielding,  pass  through  the  necessary  angle  in  the  detector  stack,  and  be  counted.  This  contamination 
was  found  to  be  significant  enough  ror  very  hard  spectra,  such  as  occurs  for  L  shells  between  1.1  and 
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1.7,  to  require  correction.  We  devised  a  correction  scheme  based  on  the  assumption  of  an  empty 
atmospheric  loss  cone.  The  correction  scheme  was  consistent  with  the  contamination  predicted  by  a 
computer  model  of  PROTEL  combined  with  a  Monte  Carlo  ray  tracing  code  developed  by  Hein  (1993). 
For  details  of  the  correction  scheme  see  Gussenhoven  et  al.  (1993). 

The  loss  cone  method  of  correcting  for  high  energy  proton  contamination  was  excessive  for  the  two 
lowest  channels  of  each  detector  head,  for  which  there  was  no  or  only  single  coincidence  in  the  counting 
scheme.  Therefore,  in  the  region  of  L  «  1.1  to  1.7,  a  cubic  spline  extrapolation  from  (interpolation 
between)  neighboring  channels  was  made  to  obtain  values  for  the  1 .5  and  2. 1  MeV  channels  (5.7  and  8.5 
MeV  channels).  These  corrections  were  made  on  the  model  average  values,  not  point  by  point. 

One  channel  (15.2  MeV)  had  significantly  lower  sensitivity  than  the  other  channels  and  one  channel 
was  an  overlap  channel  at  8.5  MeV.  To  avoid  both  redundancy  and  the  appearance  of  discontinuities  in 
differential  and  integral  spectra  at  certain  L  values,  the  8.5  and  15.2  MeV  channels  were  not  used. 

For  the  CRRES  proton  models,  a  data  base  of  differential  number  flux  values  on  the  magnetic 
equatorial  plane  was  created  using  each  PROTEL  data  point.  The  in  situ  values  were  mapped  point  by 
point  to  the  magnetic  equator  conserving  the  first  adiabatic  invariant  in  the  combined  1GRF8S  and  Olson- 
Pfitzer  quiet  magnetic  field  model.  The  equatorial  data  were  averaged  by  channel  for  each  leg  of  each 
orbit  in  L  shell  bins  of  extent  l/20th  Re  and  pitch  angle  bins  of  width  5a.  The  loss  cone  corrections  for 
high  energy  contamination  were  made  on  each  orbit  leg  average.  The  equatorial  data  from  all  orbits 
occurring  before  (after)  the  March  1991  storm  were  then  combined  in  the  same  L  and  pitch  angle  bins 
to  create  the  CRRES  quiet  (active)  proton  model. 

The  proton  inner  radiation  belt  at  a  given  energy  is  expected  to  fall  off  to  an  insignificant  value  at 
some  L.  However,  the  sporadically  occurring  solar  protons  take  over  the  flux  average  at  some  minimum 
L  that  represents  their  deepest  penetration  into  the  magnetosphere.  Although  we  retained  the  solar 
particle  effects  in  the  dose  models,  they  are  omitted  from  the  proton  belt  models.  We  feel  that  for 
predictive  purposes  the  solar  protons  are  better  represented  by  probability  rather  than  average  models  and 
will  be  handled  separately  at  a  later  time.  In  examining  the  array  of  average  flux  values  in  L  and  pitch 
angle  for  each  energy  channel,  it  is  generally  clear  where  the  fluxes  begin  to  level  off  in  L  and  transition 
from  trapped  pitch  angle  distributions  to  solar  particle  event  (SPE)  isotropic  distributions.  Identification 
of  the  boundary  between  the  two  populations  was  done  "by  hand"  for  each  energy  channel  of  each  model. 
All  but  the  1.5  MeV  channel  had  cutoff  L  values  below  5.50  R^  Thus,  5.50  is  the  maximum  L  value 
of  the  CRRES  inner  belt  proton  models. 

To  compute  the  total  proton  fluence  for  an  arbitrary  satellite  orbit,  it  is  necessary  to  completely  fill 
the  B/Bu  and  L  space  out  to  L  ■  5.50  with  omnidirectional  fluxes  at  e.'ch  grid  point  using  the  equatorial 
model  data.  We  retain  the  bln  width  of  1  /20th  R,.;  in  L  and  pioject  down  the  field  line  to  B/B„  values 
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corresponding  to  2*  steps  in  magnetic  latitude  in  a  dipole  magnetic  field  from  0°to  68°.  This  leads  to 
an  array  size  of  90  (L  shell)  by  34  (B/Bj). 

In  summary,  we  have  created  two  proton  models  from  the  CRRBS  PROTEL  data  base:  a  quiet 
model  using  data  acquired  before  the  March  1991  storm  and  an  active  model  using  data  acquired  after 
the  storm.  Each  model  potentially  consists  of  twenty-four  sub-models,  one  for  each  energy  channel.  As 
discussed  above,  we  eliminate  one  of  the  8.5  MeV  channels  and  the  15.2  MeV  channel.  This  results  in 
44  different  flux  model  files,  which  are  stored  in  binary  format  in  "FLUX_MOD.OO"  through 
"FLUX_MOD.43".  A  45th  model  file  representing  the  NASA  AP8MAX  integral  flux  for  protons  with 
energy  >90.4  MeV  is  also  stored  in  the  same  format  and  called  "NASA_MOD.OO".  This  file  is  used 
to  determine  the  integral  omnidirectional  fluence  for  the  protons  with  energies  over  90.4  MeV,  which 
is  the  upper  limit  of  the  integral  fluence  calculations  from  the  CRRES  models.  The  integration 
parameters  are  discussed  in  Section  3,  "Methods  of  Calculation". 


3.  METHODS  OF  CALCULATION 


PROTEL’s  24  energy  channels  have  the  following  lower  and  upper  boundaries  and  midpoint  energies. 
The  flux  models  are  stored  In  the  FLUXJMOD.M  files  indicated  in  the  two  columns  on  the  right. 

Energy  Channel  Parameters  (MeV)  FLUX_MOD M  file 


Channel 

E,0 

E.nd 

E„, 

Quiet 

Active 

1 

1.0 

1.5 

2.1 

FLUX_MOD.OO 

FLUX_MOD.01 

2 

2.0 

2.1 

2.3 

FLUX_MOD.02 

FLUX_MOD.03 

3 

2.2 

2.5 

2.8 

FLUX_MOD,04 

FLUX_MOD.05 

4 

2.8 

2.9 

3.2 

FLUX_MOD.06 

FLUX_MOD.07 

5 

3.1 

3.6 

4.1 

FLUX_MOD.08 

FLUX_MOD,09 

6 

3.9 

4.3 

4.8 

FLUX_MOD,10 

FLUX_MOD.ll 

7 

4.6 

5.7 

7.0 

FLUX_MOD.12 

FLUX_MOD.13 

8 

7.3 

8.4 

9.4 

no  FLUXJV10D,##  files 

9 

6.0 

6.8 

7,7 

FLUX_MOD.14 

FLUX_MOD.15 

10 

7.5 

8.5 

9.6 

FLUX_MOD.I6 

FLUX_MOD.  17 

11 

9.3 

9.7 

10.2 

FLUXJVIOD.18 

FLUX_MOD.  19 

12 

9.9 

10.7 

11.5 

FLUX_MOD.20 

FLUX_MOD.21 

13 

11.2 

13.2 

15.2 

FLUX_MOD.22 

FLUXMOD.23 
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Energy  Channel  Parameter*  (MeV) 
Channel  E* 


FLUXMOD.##  file 
Quiet  Active 

no  FLUX  MOD.##  files 


14 

14.7 

15 

15.5 

16 

18.0 

17 

25.3 

18 

26.1 

19 

34.9 

20 

37.8 

21 

44.3 

22 

53.3 

23 

62.1 

24 

73.1 

Em 

15.2 

15.9 

16.9 

18.3 

19.4 

20.8 

26.3 

27.2 

30.9 

35.6 

36.3 

37.7 

41.1 

48.1 

47.0 

53.5 

55.0 

62.1 

65.7 

73.1 

81.3 

100.0 

FLUX.MOD.24 
FLUX_MOD.26 
FLUX_MOD.28 
FLUX_MOD.30 
FLUX.MOD.32 
FLUXMOD.34 
FLUX_MOD.36 
FLUX.MOD.38 
FLUX_MOD.40 
FLUX  MOD.42 


FLUX_MOD.25 
FLUX_MOD.27 
FLUX.MOD.29 
FLUX_MOD.31 
FLUX.MOD.33 
Fl.UX_MOD.35 
FLUX.MOD.37 
FLUX.MOD.39 
FLUX_M0D.41 
FLUX  MOD.43 


Channels  8  and  14  were  omitted  from  r  !1  calculations. 

The  models  give  differential  omnidirectional  flux.  To  calculate  differential  omnidirectional  fluence 
per  year  for  an  orbit  input  by  the  user,  the  time  in  seconds  spent  in  each  (L,  B/Bb)  bin  is  calculated  by 
7IMEB1NS.EXE,  and  the  differential  flux  for  each  bin  is  multiplied  by  the  time  in  seconds  spent  in  that 
bin.  The  figures  for  the  individual  bins  are  summed.  The  resulting  figure  is  then  multiplied  by  (seconds 
per  year  /  sum  of  time  in  seconds  for  all  bins). 

Integral  omnidirectional  fluence  is  J isJ{E)dE,  where  J(E)  is  the  differential  omnidirectional 
fluence  for  energy  E.  The  integration  is  approximated  using  rectangular  steps,  so  that  integral 

eh20 

omnidirectional  fluence  *  JT  FiE^  LEi  ,  where  F(E,)  is  the  differential  omnidirectional  fluence  at 

the  average  energy  of  channel  ch,.  Because  of  overlapping  energy  ranges,  channels  5  and  15  were 
omitted  from  the  integral  omnidirectional  fluence  calculation  (in  addition  to  channels  8  and  14  which  are 
omitted  from  all  calculations)  and  new  boundaries  were  set  up  to  prevent  multiple  contribution*  from  the 
same  energy.  The  integration  summation  for  a  channel  ch,  begins  at  the  lower  boundary  of  that  channel 
and  ends  at  the  new  upper  boundary  of  the  highest  energy  channel,  which  is  now  channel  20  with  an 
upper  boundary  of  90.4  MeV.  It  is  left  for  the  user  to  interpolate  between  energies  if  desired.  The 
energy  bounds  and  AEs  used  in  the  integrations  follow. 
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Integration  Parameters  (MeV) 

Channel  # 

E* 

EmU 

E* 

AE 

1 

1.1 

1.5 

1.9 

0.8 

2 

1.9 

2.1 

2.3 

0.4 

3 

2.3 

2.5 

2,7 

0.4 

4 

2.7 

2.9 

3.1 

0.4 

5 

3.1 

4.3 

5.5 

2.4 

6 

5.5 

5.7 

5.9 

0.4 

7 

5.9 

6.8 

7.7 

1.8 

8 

7.7 

8.5 

9.3 

1.6 

9 

9.3 

9.7 

10.1 

0.8 

10 

10.1 

10.7 

11.3 

1.2 

11 

11.3 

13.2 

15.1 

3.8 

12 

15.1 

19.4 

23.7 

8.6 

13 

23.7 

26.3 

28.9 

5.2 

14 

28.9 

30.9 

32.9 

4.0 

15 

32.9 

36.3 

40.2 

7.3 

16 

40.2 

41.1 

43.2 

3.0 

17 

43.2 

47.0 

30.8 

7.6 

18 

50.8 

55.0 

59.2 

8.4 

19 

59.2 

65,7 

72.2 

13.0 

20 

72.2 

81.3 

90.4 

18.2 

4.  LOGICAL  FLOW  OF  CRRESPRO 


Figure  2  shows  a  flow  diagram  of  the  CRRESPRO  execution.  The  solid  arrows  indicate  the  direction 
of  flow,  and  the  dashed  arrows  indicate  the  input  and  output  at  each  step.  After  each  logical  step, 
CRRESPRO  stops  and  presents  the  user  with  information  about  what  it  will  be  doing  next,  then  waits  for 
a  key  to  be  pressed  before  continuing. 

First  the  user  must  enter  the  name  of  the  orbit.  This  name  will  be  used  with  all  the  flies  generated 
by  CRRESPRO.  No  extension  can  be  given  with  the  name,  but  a  path  can  be  specified  with  the  name 
if  the  flies  are  to  be  written  to  some  other  directory.  In  Figure  2,  the  input  name  is  represented  by 
"7NAME?"  First  time  users  of  CRRESPRO  should  type  "CRRESPRO"  on  the  command  line.  CRRESPRO 
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rXRX  CRRESF 

~z 

Input  n*m«  of 
run: 

Input  orbital 


Umt  input 
Uaer  input 


- - * 

don  van  to  moan  I 

orfal  tn 

IjlammaJ 

—  Umt  input 
i  ^Tnnmrf.ELX 


—  7mo»?.EUC 


1  action? 


DtPOU.EXE 


CRRESF 


—  7MOW7.SPH 


—  7uum?.MAO 

TIM 

FLUX_MOD.fi 

—  Tnamrt.TIM,  7nam«?.ELf 
■  '►7nnnu7.DAT 


—  Tnamtf.DAT 

.  —  7msm7.DAT 

FLUX.MOD.li 
_  7iuum?.MA0 


7ruaw?.SPH 
’  7ruum7.DAT 


Figure  2.  Simplified  Flow  Diagram  of  CRRESPRO.  Solid  Arrows  Refer  to  Logical  Flow,  and  Dashed 
Arrows  Refer  to  Input  and  Output  at  Each  Step. 


will  query  the  liter  for  the  name  when  it  runt.  (If  filet  from  a  previous  CRRESPRO  run  with  the  tame 
name  exitt,  CRRESPRO  will  ask  the  user  if  these  files  should  be  used.)  A  second  option  for  users 
familiar  with  CRRESPRO  is  to  include  the  name  on  the  command  line  by  typing  "CRRESPRO 
?NAME?.n  When  this  option  is  used,  CRRESPRO  does  not  stop  to  tell  the  user  what  it  is  doing  or  wait 
for  user  input.  This  allows  long  runs  taking  many  hours  to  be  done  without  requiring  the  user  to  respond 
to  the  CRRESPRO  prompts.  In  this  mode,  any  intermediate  files  that  exist  with  the  given  name  will 
automatically  be  used.  Note:  If  a  previous  CRRESPRO  orbit  of  the  same  name  was  interrupted  during 
execution  and  the  resulting  files  still  exist,  an  error  will  result.  The  user  should  remove  the  old  files  or 
choose  a  new  name  for  the  orbit. 

Next,  CRRESPRO  prompts  the  user  for  the  orbital  elements.  These  elements  describe  the  orbit  to 
be  tested.  Only  orbits  for  times  from  day  1 , 1985  to  day  365, 1999  can  be  used.  There  are  three  options 
for  specifying  orbital  elements: 

Position  and  Velocity;  The  elements  for  position,  f„  ry,  r„  and  velocity,  v,(  vy,  v„  are 
specified  in  km  and  km/sec  using  the  Earth  Centered  Inertial  (ECI)  coordinate  system.  Time 
is  given  in  year,  day  of  year,  hours,  minutes,  and  seconds  and  must  be  fall  between  the  start 
of  day  l,  1985  and  the  end  of  day  365,  1999. 

Mean  Orbital  Elements:  Inclination,  right  ascension  of  the  ascending  node,  argument  of 
perigee,  and  mean  anomaly  are  specified  in  degrees,  mean  motion  is  specified  in  revolutions 
per  day,  and  eccentricity  is  specified.  Time  Is  given  in  the  same  format  as  for  the  position  and 
velocity  elements. 

Solar  Elements:  Inclination  is  specified  in  degrees,  while  altitudes  of  apogee  and  perigee  are 
specified  in  km  assuming  Re  -  6378  km,  and  local  times  of  apogee  and  maximum  inclination 
are  specified  in  hours.  These  elements  are  useful  to  scientists  and  engineers  concerned  with 
the  location  of  an  orbit  with  respect  to  the  sun  and  to  those  who  are  not  familiar  with  the 
standard  orbital  elements.  The  solar  elements  are  not  a  full  set  of  elements  because  there  are 
only  five  elements  (besides  time)  to  input.  The  sixth  element  results  from  the  assumption  that 
the  orbit  is  at  apogee  at  the  time  for  which  the  elements  are  specified.  Local  time  of  apogee 
of  12.0  indicates  that  the  satellite  is  at  apogee  directly  between  the  sun  and  Earth.  Local  time 
of  maximum  inclination  indicates  the  local  time  where  the  geographic  latitude  of  the  satellite 
is  at  a  maximum  positive  value  Default  values  of  12.0  and  18.0  are  used  respectively  for  the 
two  local  times,  but  these  values  may  be  changed.  A  six  hour  difference  between  the  local 
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time  of  apogee  and  the  local  time  of  maximum  inclination  causes  apogee  to  be  at  0*  geographic 
latitude.  Oftentimes  this  arrangement  results  in  the  highest  fluence  since  the  highest  proton 
concentration  is  near  the  magnetic  equator,  close  to  the  geographic  equator.  There  are  other 
cases  where  the  6-hour  difference  does  not  result  in  the  highest  fluence.  One  example  is  an 
elliptical  polar  orbit  with  apogee  in  the  equatorial  slot  region.  The  format  of  the  given  time 
is  as  specified  for  the  position  and  velocity  elements.  The  solar  elements  cannot  be  used  by 
LOKANGL.EXE  directly,  so  CRRESPRO  converts  solar  elements  into  mean  orbital  elements. 
(LOKANGL.EXE  is  described  below;  see  Appendix  C  for  calculations.) 

The  elements  are  checked  to  make  sure  the  orbit  is  not  hyperbolic  and  perigee  is  not  below  200  km  (see 
Appendix  C).  The  time  step  used  in  the  study  is  chosen  so  that  the  orbit  will  sweep  -2°  arc  at  perigee 
in  a  single  time  step  (see  Appendix  C).  The  user  then  inputs  the  start  time  (between  day  1  1985  and  day 
365  1999)  and  length  of  the  study.  The  elements  are  written  to  an  ASCII  file  called  7NAME7.ELE. 
After  the  user  is  done  entering  the  set  of  elements,  CRRESPRO  waits  for  HD"  to  be  pressed. 

CRRESPRO  calls  LOKANGL.EXE,  which  reads  the  element  file  and  generates  an  ephemeris  file, 
7NAME7.SPH,  which  contains  geographic  position.  The  position  is  calculated  at  each  time  increment 
in  the  specified  period  of  study.  Appendix  D  contains  a  description  of  the  7NAME7.SPH  file. 
LOKANGL.EXE  is  used  at  the  Phillips  Laboratory  Geophysics  Directorate  to  predict  positions  of  Barth 
satellites  using  orbital  elements  as  an  input.  The  code  was  modified  by  RADEX  Inc.  for  use  on  the  IBM 
PC.  Contact  this  office  (PL/GPSP)  for  information  on  the  format  of  the  7NAME7.ELE  file  as  well  as 
a  detailed  description  of  LOKANGL.EXE.  LOKANGL.EXE  calculates  the  position  of  an  Earth  satellite 
using  a  perturbation  expansion  with  the  second  and  third  zonal  harmonics  of  Earth’s  gravitational  field. 
It  is  not  considered  accurate  for  times  more  than  100  days  from  a  specified  orbital  element. 
LOKANGL.EXE  cannot  evaluate  hyperbolic  or  parabolic  orbits  or  orbits  where  the  radius  becomes  less 
that  one  Earth  radius. 

When  LOKANGL.EXE  is  finished,  CRRESPRO  will  give  the  user  the  option  to  view  the  orbit  by 
calling  ORBIT.EXE.  Appendix  E  contains  a  description  of  ORBIT.EXE.  If  the  file  name  was  given  on 
the  command  line,  CRRESPRO  will  not  give  the  option  to  view  the  orbit  and  will  proceed  on  to  the  next 
step. 

CRRESPRO  calls  MAGMODEL.EXE  (see  Appendix  F)  to  convert  the  geographic  coordinates  in 
7NAME?  SPH  into  magnctospheric  coordinates  that  are  put  into  the  binary  file  7NAME7.MAG,  also 
described  in  Appendix  D.  The  routines  used  by  MAGMODEL.EXE  to  calculate  L  shell  and  B/Bt,  were 
developed  for  the  Air  Force  by  Karl  Pfitzer  at  McDonnell  Douglas  Space  Systems  Corp  (PfUzer,  1991). 
The  IGRF85  model  used  by  MAGMODEL.EXF  is  only  valid  for  times  after  the  start  of  day  1,  1985. 
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TIMEBINS.EXE  is  called  next  by  CRRESPRO.  This  program  reads  the  data  in  the  magnetospheric 
ephemeris  file  and  interpolates  between  points  using  a  cubic  spline  to  determine  how  much  time  the  orbit 
spends  in  each  of  the  bins  defined  for  the  flux  models.  The  resulting  array  is  written  to  the  binary  file 
?NAME?.TIM. 

CRRESPRO  then  calls  PLUXDATA.EXE.  This  program  reads  in  the  times  in  each  bin  from 
7NAME7.TIM  and  calculates  the  differential  and  integral  omnidirectional  ftuence  per  year  as  explained 
in  Section  3,  "Methods  of  Calculation".  PLUXDATA  assumes  that  the  first  44  files  (##  »  00-43)  are 
the  standard  flux  files;  these  flies  should  not  be  renamed.  Any  other  FLUX_MODJ#  flies  found  in  the 
active  directory  will  also  be  integrated  and  included  in  the  output.  The  integrated  data  is  output  to  a 
formatted  ASCII  file  called  7NAME7.DAT.  Appendix  F  shows  the  contents  of  SAMPLE.DAT  which 
was  generated  in  Appendix  A.  The  top  section  gives  the  name,  orbital  elements,  and  the  times  of  the 
study.  After  the  time  of  study,  the  amount  of  the  orbit  (by  time)  not  covered  by  the  flux  bins  is  given. 
An  orbit  similar  to  the  CRRES  orbit  is  used  in  SAMPLE.DAT.  Since  the  proton  models  extend  only  to 
L  -  5.5,  part  of  the  orbit  is  out  of  range  of  the  model.  For  SAMPLE,  this  is  44  percent  (by  time),  as 
indicated.  With  the  exception  of  the  1 .5  MeV  channel,  the  flux  bins  not  covered  by  CRRESPRO  would 
make  no  contribution  to  the  fluence.  (See  Section  2,  "Flux  Models".)  The  next  part  of  the  CRRESPRO 
output  is  a  table  showing  differential  omnidirectional  fluence  per  year  for  the  quiet  and  active  periods. 
This  is  followed  by  a  table  showing  integral  omnidirectional  fluence  per  year  for  the  quiet  and  active 
periods.  The  first  44  FLUXMOD.M  files  are  used  to  produce  the  tables.  The  file  NASA_MOD.OO 
is  accessed  to  produce  the  final  line  of  the  CRRESPRO  output,  which  shows  the  integral  omnidirectional 
fluence  for  protons  with  energies  >90.4  MeV. 

At  this  point  CRRESPRO  is  done  with  calculations  and  provides  several  options  to  view  the  data. 
Pressing  "P"  calls  the  DOS  command  "COPY  7NAME7.DAT  PRN,"  which  will  print  the  ASCII  file 
7NAME7.DAT.  Pressing  "V"  displays  the  file  on  the  screen.  "G"  graphs  the  orbit  using  ORBIT.EXE. 
"D"  graphs  the  orbit  superimposed  on  the  flux  models  using  DIPOLE.EXE  (see  Appendix  E).  and  "Q" 
quits  CRRESPRO.  A  benchmark  file  is  created  by  CRRESPRO  which  shows  how  long  each  of  the 
different  utilities  took  to  run.  This  file  has  the  name  7NAME7.BEN  and  is  discussed  in  more  detail  in 
Appendix  A. 
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Appendix  A 


Sample  Session 


Type: 

CRRESPRO  <  Enter  > 

Press  any  key  once  you  have  read  the  information  on  the  screen.  You  are  now  asked  for  a  file 
name.  Here  we  use  "Sample"  as  the  file  name.  Type: 

Sample  <  Enter  > 

Press  any  key  once  you  have  read  the  information  on  the  screen.  Press  "S"  to  select  solar  elements. 
On  the  next  screen  we  specify  an  orbit  that  is  very  similar  to  the  orbits  used  by  CRRES.  Enter  the  solar 
elements  and  the  time  for  which  those  elements  apply: 

year  -  1991 
day  -  130 
hour  ■  8 
minute  «  0 
second  *  0 
Inclination  «  18 
Perigee  ■»  350 
Apogee  “  33000 
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Leave  the  two  default  local  times  as  they  are  at  12  and  18.  This  causes  the  satellite  to  be  at  apogee 
at  12:00  and  at  maximum  inclination  at  18:00  local  time.  (See  Section  3  for  the  applicable  logic.)  If  you 
enter  something  incorrectly,  you  can  tab  back  to  that  line  and  enter  it  again.  Pressing  "x"  then  <  Enter  > 
on  a  line  will  delete  it.  When  you  are  done  press  "D."  At  the  next  screen  enter  the  time  of  the  study: 

year  »  1991 
day  -  130 
hour  »  8 
minute  -  0 
second  «  0 
length  -  1  (days) 

The  time  of  the  study  can  differ  from  the  time  entered  previously  for  the  element  information.  Press 
"D"  when  you  are  done. 

A  message  will  appear  on  the  screen  estimating  the  calculation  time  for  a  66  Mhz  80486  computer. 
In  this  case  the  fluence  will  take  about  1.1  minutes  to  calculate.  The  actual  time  required  will  vary  with 
the  speed  of  the  computer  and  the  type  of  orbit.  In  general,  lower  altitudes  take  longer  than  higher 
altitudes,  elliptical  orbits  take  longer  than  circular,  and  high  inclination  orbits  take  longer  than  low 
inclination.  If  you  wish  to  terminate  the  program  at  this  point  press  "N;"  otherwise  press  "Y"  to 
continue,  and  then  press  any  key  once  you  have  read  the  information  on  the  screen. 

The  code  LOKANOL.EXE  will  execute.  In  the  upper  left  corner  of  the  screen  is  written  "00001 . 1  ” 
which  is  the  estimated  time  to  completion  in  minutes.  Below  it  is  written  the  actual  execution  time.  Once 
LOKANGL.EXE  is  done  executing,  CRRESPRO  will  ask  if  you  would  like  to  see  the  orbit.  Press  "Y” 
to  execute  ORBIT.EXE  which  shows  a  view  of  the  orbit.  The  view  in  the  lower  right  looks  down  on 
Earth’s  north  pole  with  the  sun  towards  the  top.  The  other  two  images  show  side  views,  looking  along 
X  and  Y.  The  orientation  of  the  axes  on  these  views  follows  standard  drafting  conventions.  Press  “Q1' 
to  quit  ORBIT.EXE,  then  press  any  key  after  you  have  read  the  information  on  the  CRRESPRO  screen. 

The  code  MAOMODEL.EXE  will  execute.  The  actual  execution  time  (upper  left  corner  of  screen) 
does  not  start  from  zero,  but  begins  where  it  left  off  when  LOKANGL.EXE  was  finished  executing. 
When  MAGMODEL.EXE  is  through,  press  any  key  after  reading  the  information  on  the  screen. 

The  code  TIMEBINS.EXE  will  execute.  This  code  runs  faster  than  LOKANGL.EXE  or 
MAGMODEL.EXE.  When  it  is  done,  press  any  key  after  reading  the  information  on  the  screen. 

The  code  FLUXDATA.EXE  will  execute.  The  time  taken  by  this  code  does  not  depend  on  the  size 
of  the  run.  When  it  is  done,  read  the  information  on  the  screen  then  press  any  key. 
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You  now  have  several  options  for  viewing  the  data.  If  you  have  a  printer,  press  "P"  to  get  a 
printout.  (Appendix  F  is  a  copy  of  the  printout  of  the  file  SAMPLE.DAT.)  You  can  press  "V"  to  see 
the  output  on  the  screen.  "Q"  will  return  you  to  the  menu.  Fressing  "G"  calls  ORBIT.EXE  to  show  the 
orbit  plot  again.  Notice  that  more  information  is  displayed  in  the  upper  left  corner  now  that  the  job  is 
complete.  This  information  comes  from  the  header  in  SAMPLE.DAT.  "D"  executes  DIPOLE.EXE 
which  shows  the  orbit  track  through  the  radiation  belts.  The  orbit  track  is  plotted  entirely  in  the  upper 
quadrant.  Although  this  may  at  first  seem  unnecessarily  confusing,  it  does  allow  one  to  use  the  lower 
quadrant  to  view  the  radiation  belt  intensities  without  orbit  track  obstruction.  Notice  that  for  "Sample. " 
part  of  the  orbit  is  outside  the  high  L  limit  of  5.5.  With  the  exception  of  the  1.5  MeV  channel,  this  does 
not  mean  that  the  calculation  is  missing  protons  from  the  radiation  belts.  The  statement  on  the  printout 
indicating  that  44.0  percent  of  the  orbit  is  out  of  range  reflects  the  time  spent  outside  the  maximum  L 
•  5.5,  maximum  magnetic  latitude  -  68°  boundaries  of  the  model. 

The  time  required  for  your  computer  to  predict  the  fluence  for  orbit  "Sample"  is  recorded  in  the 
ASCII  benchmark  file  called  SAMPLE. BEN.  Following  is  a  copy  of  the  bench  mark  file  for  Sample 
used  on  a  66  Mhz  486  GATEWAY  computer  at  the  Phillips  Laboratory. 

*****  Time  increment  -  136.34  seconds. 

13:47:53.40 ->  LOKANGL  started 
13:47:55.98  -->  LOKANGL  finished 
Number  of  points  -  635 
13:47:56.09  ~>  MAGMODEL  started 
13:48  04.28  -•>  MAGMODEL  finished 
13:48:04.28  -->  TIMEBINS  started 
13:48:04.66  ~>  TIMEBINS  finished 
13:48:04.72  ~>  FLUXDATA  started 
13:48:24.82  »>  FLUXDATA  finished 

Time  increment  is  the  step  size  chosen  by  CRRESPRO  to  calculate  ephemeris  and  fluence.  This 
time  was  chosen  so  that  the  satellite  swings  through  approximately  2e  arc  at  perigee.  See  Appendix  C 
for  calculations.  LOKANGL  took  2.58  seconds  to  run,  MAGMODEL  took  8.19  seconds,  TIMEBINS 
took  0.38  second,  and  FLUXDATA  took  20. 1  seconds.  The  rest  of  the  time  was  spent  waiting  for  user 
input  and  was  not  considered  in  the  estimated  run  time  of  1. 1  minutes.  The  actual  time  taken  was  30.49 
seconds  while  the  estimate  was  for  1.1  minutes.  The  estimate  tends  to  be  a  worst  case  estimate.  By 
comparing  your  benchmark  file  to  the  one  shown  above  you  can  calculate  the  speed  of  your  computer 
versus  the  benchmark  computer.  This  ratio  applies  to  the  estimated  run  time. 
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Appendix  B 

Format  of  FLUXJV10D.##  Files 

The  FLUX_MOD  .00  files  contain  data  binned  by  L  shell  and  B/B0.  Any  rate  data  such  as  flux  or  dose 
rate  can  be  put  into  these  tiles,  provided  the  following  format  is  used: 

-  The  filename  must  be  FLUX_MOD.##  (00  -  00  to  99) 

-  The  file  structure  is  of  the  type: 

-  BORLAND  PASCAL 


FluxRecotd  -  record 
MainTitle, 

SubTitle, 

Units  :  String[136); 

d  :  Array[1..90,l.. 34]  of  Single; 


-  ANY  OTHER  LANGUAGE 

MainTitle, 

SubTitle, 

19 


Units 


Each  is  an  array  of  136  eight  bit  bytes  or  characters.  The  first  byte 
contains  the  number  of  characters  in  the  array  (0  to  136).  The  following 
bytes  contain  the  characters  in  the  array, 
d  This  is  the  data.  Each  bin  is  a  4  byte  real  number  following  the  convention 

set  by  the  Intel  math  coprocessor.  The  data  is  binned  in  a  two  dimensional 
array  of  1  to  90  (L  shell)  by  1  to  34  (B/Bq).  Consecutive  B/B„  bins  occupy 
consecutive  spaces  in  memory. 

•  The  files  have  four  different  groups  of  information. 

-  The  first  three  groups  are  strings.  MainTitle  and  SubTitle  are  strings  of  characters  that  are 
displayed  as  the  two  title  lines  by  the  utility  D1P0LE.EXE.  Units  is  the  string  displayed  by 
DIPOLE.EXE  near  the  scale  to  identify  what  type  of  data  is  displayed. 

-  The  fourth  group  is  the  data,  d  is  a  two  dimensional  array  of  1  to  90  (L  shell)  by  1  to  34 
(B/B0)  bins. 

--  i  is  the  index  associated  with  L  shell,  and  ranges  from  1  to  90.  The  limits  of  bins  with 
index  I  are:  0.95  +0.051  RE  £  L  shell  <  1.00  +0.051  RE 

—  J  is  the  index  associated  with  B/B0,  and  ranges  from  1  to  34.  The  limits  of  bins  with 
index  j  are: 

Bh  S  B/B0  <  Bjj 

B0  j.,  -  1.000,  1.004,  1.020,  1.046,  1.085,  1.140,  1.200,  1.300,  1.400,  1.520, 
1 .690,  1 .880, 2. 100,  2.400,  2.730,  3.130,  3.670, 4.350,  5.020,  6. 100,  7.410,  9.088, 
1 . 1 29e 1 ,  1.422el,  1.816el,  2.356el,  3.107el,  4.147el,  5.723el,  8.025el,  1.154e2 
1.707e2,  2.607e2,  4.134e2  ,  6.846e2.  (where  e#  indicates  (  10^).) 

—  The  values  of  B,  are  such  that  the  limits  B0  to  B34  c«.  noroximately  ±68°  magnetic 
latitude  <n  a  dipole  field. 
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Appendix  C 

Calculations  Used  to  Create  Element  File 


We  give  here  the  equations  used  to  1)  cp'culatc  the  time  step  used  in  LOKANGL.EXE,  2)  convert  solar 
elements  to  mean  orbital  elements,  and  3)  eliminate  invalid  orbits,  namely  those  that  are  not  closed  or 
that  pass  close  to  or  beneath  the  surface  of  Earth.  Two  constants  specific  to  Earth  that  are  used 
throughout  are  /*,  tKe  gravitational  constant  times  the  mass  of  Earth,  and  RE,  the  radius  of  Earth  assuming 
Earth  is  spherical.  These  values  are: 

H  -  398601.2  km’sec1, 

Re  -  6378.145  km. 

U)  Tims  Step  Calculation; 

If  Position/Velocity  Elements  are  used  where 

r  -  \jr\  +  r*  +  /•*  ,  v  .  ftjTvf+v2', 

then  define  the  three  quantities,  a,  h,  and  e  as 

a-'.l-*  .  h  .  |7xT|  ,  .  1- *1  . 

r  u  na 

Here  a  is  the  semimajor  axis  in  km,  h  is  the  angular  momentum  in  kmVsec,  and  e  is  the  eccentricity. 
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If  Mean  Orbital  Elements  are  used,  the  eccentricity  t  and  the  mean  motion  Mm  in  revolutions  per 
day  are  given  as  elements.  The  semimajor  axis  a  and  the  magnitude  of  the  angular  momentum  h  are 
given  by: 

o5  -  mK12‘3600)/(»A/JJ1  ,  h1  «  fia(l  -eJ) 

In  terms  of  a ,  h,  and  e,  the  time  step  T„  which  is  approximately  equal  to  the  time  to  traverse  2°  arc  at 
perigee  is: 

T,  -  ai(l-es)/r'^.  . 

In  this  calculation  d@/dt  at  perigee  is  assumed  to  be  constant  over  the  2 9  arc. 


Use  the  convention  day  «■  1  on  1  January  1992. 

Define  the  following  quantities: 
r,  ■  A,  +  Rb; 

rf  -  Ap  +  Re;  A,  and  Ap  are  altitudes  from  the  surface  of  Earth  in  km; 

He  -  6.594703  +0.0657098243^+ 1.00273791 'Hr  [U.S.  Naval  Observatory.  1992J; 
0,  «=  (H,  +  12)*W12,  where  H,  is  local  time  of  apogee; 

0,  =  (H,  -  6)*r/12,  where  H,  is  local  time  of  maximum  inclination. 

Then  the  six  elements  are: 

inclination,  i,  is  given; 
eccentricity,  e  »  (r,  -  rp)  /  (r,  +  r^; 

mean  anomaly,  M  =  180°  (satellite  is  at  apogee  when  elements  are  defined); 

squared  mean  motion.  na  ■  n/[ (r, + rp)]3(  1 2  •3600/t)2; 

right  angle  of  ascending  node,  0  *  15(H,  -  6  +  He); 

argument  of  perigee,  w,  cos(w)  »  cos(0t)cos(0,)  +  sin(0,)sin(0,). 


The  orbit  is  valid  if: 


a)  the  orbit  is  closed:  0  s  e2  <  1; 

b)  perigee  is  not  below  the  surface  of  Earth:  e(l-e)  *  6400  km. 
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Appendix  D 


Format  of  Ephemeris  Files 


The  geographic  ephemeris  files  end  with  the  extension  "SPH."  These  are  binary  files  that  are  output 
from  LOKANGL.EXE.  Each  record  of  the  file  specifies  the  position  of  a  satellite,  and  the  time  interval 
between  all  records  is  constant  (T,  is  in  Appendix  C).  These  files  are  input  by  MAGMODEL.EXE  to 
produce  the  magnetospheric  ephemerides. 

Each  record  of  the  file  contains,  in  the  following  order: 


year 

2  byte  integer 

day 

2  byte  integer 

second 

4  byte  real 

latitude 

4  byte  real  (deg) 

longitude 

4  byte  real  (deg) 

radius 

4  byte  real  (km) 

The  magnetospheric  ephemeris  files  end  with  the  extension  "MAG."  These  are  binary  files  that  are 
output  from  MAGMODEL.EXE.  Each  record  of  the  file  specifies  the  position  of  a  satellite  in  Earth’s 
magnetic  field,  and  the  time  interval  between  the  records  is  constant  (T,).  The  requirement  for  a  constant 
time  step  is  set  by  TIMEB1NS.EXE.  These  files  are  input  by  TIMEBINS.EXE  to  produce  files  indicating 
how  much  time  the  orbit  spends  in  the  magnetospheric  bins  defined  for  the  flux  models  (see  Appendix 
B). 

Each  record  of  a  ".MAG"  file  contains,  in  the  following  order: 

L  shell  :  4  byte  real  (Re) 

B/B„  :  4  byte  real 

Open  field  lines  are  indicated  by  an  L  shell  of  -1  or  L  shells  >  15  Rr. 
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The  last  record  of  each  file  contains,  in  the  following  order: 

•100  :  4  byte  real 

TimeStep  :  4  byte  real  (time  interval  between  steps) 

LOKANGL.EXE  can  be  bypassed  by  creating  a  ".SPH"  ephemeris  file  of  the  proper  format. 
CRRESPRO  will  still  request  orbital  elements  so  you  must  generate  a  dummy  element  file,  but 
CRRESPRO  will  not  execute  LOKANGL.EXE  if  the  ".SPH"  file  exists.  In  the  same  way,  a  different 
magnetospheric  model  can  be  used  by  creating  a  '.MAG"  file  of  the  proper  format.  CRRESPRO  will 
still  request  orbital  elements,  and  will  run  LOKANGL.EXE  (if  a  ".SPH"  file  does  not  exist),  but  it  will 
not  execute  MAGMODEL.EXE. 
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Appendix  E 

CRRESFRO  Support  Programs 


Ooordinat*  swat an  < Barth  oantarad>: 
Z  par* Hal  to  Earth'*  spin  axis 
V  parp.  to  2  «  Earth-tun  lino 
X  font*  ri*ht  handad  coordinator 


•  040.00000  kn 


dau  170  IMS  12:00:00.00 
nit  nooo  -  • *0.00000  kn 
nit  Pari 
Inal inat 
nr*  Pori 

BAA  M - 1  — 

M  nnonal 

Eoaantr 1 
H  Notion 


100.00000  da* 


100.00000  da* 

0.00000 
14.10070  rod 


Start  tin*  of  study: 

daw  170  lOtt  00:00:00.00 
Stop  tin*  of  atudu: 
day  170  lOOO  10:00:00.00 
_ Praaa  ,’QM  to  quit 


Dunk 


-V 


Sou 


North 


-V 

Dawn 


Ni*ht 


Figure  El.  ORBIT  EXE  Screen 


ORBIT.EXE  -  This  is  a  display  utility  program  that  works  on  a  color  EGA/VOA  monitor.  It  is 
called  by  CRRESPRO,  but  may  also  be  called  from  the  command  line.  This  program  reads  the  ".SPH" 
files,  and  will  also  read  the  ".DAT"  files  if  they  exist.  When  calling  ORBIT.EXE  from  the  command 
line,  the  name  of  the  run  must  be  included  on  the  command  line.  The  name  of  these  files  may  include 
a  path,  but  not  an  extension.  If  one  typed  "ORBIT  \ORB\SAMPLE,"  the  data  file  SAMPLE. SPH  in 
the  directory  "\ORBN  would  be  displayed.  If  SAMPLE.DAT  exists  then  orbital  elements  read  from  that 
file  would  be  displayed  in  the  upper  left  comer.  Figure  El  shows  a  screen  capture  of  ORBIT.EXE  using 
the  files  created  in  the  sample  session  of  Appendix  A.  The  upper  left  comer  shows  information  about 
the  orbit  drawn,  In  the  coordinate  system  used  (see  Figure  E2),  Z  points  along  Earth's  spin  axis  and  is 
positive  toward  the  north  pole,  while  Y  is  perpendicular  to  Z  and  the  Earth-Sun  line  and  is  positive 
towards  dusk.  X  completes  the  right  handed  coordinate  system  and  goes  through  the  noon  meridian.  The 
lower  right-hand  comer  of  Figure  El  shows  the  orbit  projected  into  the  XY  plane.  Tick  marks  show 
integer  numbers  of  Re.  The  other  two  comers  show  the  orbit  projected  into  the  XZ  and  YZ  planes.  The 
orientation  of  the  coordinates  follows  drafting  convention  where  the  figure  in  the  lower  right  comer  is 
the  top  view  and  the  other  two  figures  are  side  views.  Pressing  "Q"  exits  the  program. 

MAGMODEL.EXE  -  This  program  uses  the  IGRF85  and  Olson-Pfitzer  Static  Magnetic  Field  Model 
IBM  PC  Fortran  routines  developed  by  Karl  Pfitzer  (1991).  To  use  this  code  two  command  line 
parameters  are  required.  The  first  is  the  path  of  the  source  file  (including  name  and  extension),  and  the 
second  is  the  path  of  the  destination  file  (including  name  and  extension).  The  source  file  is  a  geographic 
ephemeris  file,  and  the  destination  file  is  a  magnetospheric  ephemeris  file.  The  formats  of  both  these  files 
are  described  in  Appendix  D. 

DIPOLE.EXE  -  This  program  displays  the  flux  data  in  the  FLUX_MOD.M  files  The  data  is 
placed  spatially  on  the  screen  as  it  would  look  if  Earth's  magnetic  field  were  a  dipole.  (The  image  will 
not  be  proportioned  correctly  if  .m  EGA  card  is  used).  Figure  E3  shows  a  copy  of  the  screen  when 
DIPOLE.EXE  is  displaying  the  orbit  generated  in  Appendix  A.  The  figure  is  shown  in  gray  scale,  but 
the  program  is  displayed  in  color.  The  following  description  explains  what  is  seen  here  (in  gray  scale). 
The  half  circle  on  the  left  side  of  the  screen  represents  Earth  with  the  magnetic  axis  pointing  straight  up. 
Grt/  scales  represent  intensity  of  flux.  Higher  flux  tends  toward  the  darker  shading,  and  lower  flux  tends 
toward  the  lighter  shading.  If  the  flux  is  above  the  color  scale,  it  will  be  shown  as  black,  and  if  the  flux 
is  below  the  color  scale  or  there  is  no  data,  it  will  be  shown  as  light  gray.  The  display  on  the  computer 
screen  will  show  Earth  as  a  cyan  colored  half  circle.  Color  scales  represent  the  intensity  of  flux.  Higher 
flux  tends  toward  red  and  lower  flux  tends  toward  blue.  Flux  above  the  color  scale  is  shown  in  white, 
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Figure  E2.  Coordinate  System  Used  for  ORBIT.EXE. 

and  flux  below  the  color  scale  is  shown  in  gray.  The  space  bar  calls  up  an  options  menu.  The  color 
scale  can  be  reset  and  can  also  be  changed  from  a  log  to  a  linear  scale.  One  can  move  to  a  higher 
(lower)  flux  model  by  pressing  the  plus  (minus)  key.  Pressing  "Q"  quits  the  program.  The  orbit  in  L 
shell  versus  B/B0  is  plotted  as  a  black  line  on  the  upper  quadrant  of  the  flux  data.  The  orbit  is  plotted 
in  only  one  quadrant  because  B/B0  contains  no  information  about  +/-  latitude.  This  causes  the  orbit  to 
be  folded  into  the  upper  quadrant.  In  Figure  E3,  the  orbit  covers  the  area  from  Earth’s  surface  to  beyond 
the  L  cutoff  of  5.5.  The  flux  model  in  Figure  E3  shows  omnidirectional  flux  for  26.3  MeV  protons 
during  the  active  time  period.  Two  proton  belts  with  peaks  at  1.5  and  2.3  RE  clearly  show  in  the  figure. 

The  orbit  plot  superimposed  on  the  flux  data  can  be  used  to  analyze  the  effectiveness  of  CRRESPRO 
for  the  given  orbit.  Portions  of  low  orbits  may  actually  be  drawn  below  Earth’s  surface  on  the  plots. 
This  results  from  the  fact  that  Earth  is  not  actually  spherical,  the  magnetic  field  is  not  actually  dipolar, 
and  the  center  of  the  magnetic  Held  is  not  at  the  center  of  Earth.  The  valid  range  of  the  PROTEL  data 
can  be  seen  clearly  from  looking  at  the  flux  plots. 
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OMNI-DIKBCTIONAI  FLUX  —  ACUTE 
36.3  MeV 


FUe  36  limits  o I  CKRESFRQ  Model 

Figure  E3.  DIPOLE.EXE  Screen  Displaying  SAMPLE. MAO  Orbit. 

DIPOLE.EXE  can  be  run  from  the  command  line  by  typing  "DIPOLE."  Dipole  will  ask  for  the 
file  number  to  be  displayed,  and  will  then  graph  the  data.  The  file  number  can  also  be  specified  on  the 
command  line.  To  show  FLUX_MOD,09  one  would  type  "DIPOLE  9."  A  second  command  line 
parameter  can  be  included  to  show  an  orbit  superimposed  on  the  flux  map.  The  second  parameter  must 
be  the  path  (including  name  and  extension)  of  a  magnetospheric  ephemeris  file.  The  file  must  be  of  the 
format  given  in  Appendix  D. 
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#/cm*/s/IfeV 


Appendix  F  | 

j 

Contents  of  SAMPLE.DAT 


CRRB8PR0  Prediction  Utility  >  lAMPLg 


Element  1 1  (Mi  day  130  1991  07i59i59.71 
Alt  Apog  ■  32999.99933  km  Alt  Perl  *  390.00066  km 

Inellnnt  ■  18.00000  deg  Arg  Perl  ■  180.00000  deg 

RAA  Node  ■  47.29340  deg  N  Anonal  ■  180.00000  deg 

Eceentr!  ■  0.70819  H  Notion  >  2.48032  rpd 


8tert  time  of  etudyi  dey  130  1991  O81OO1OO.OO 
stop  time  of  etudyi  dey  131  1991  OOiOOtGO.OO 
Time  of  etudy  It  0.99990  deye  <44. OX  of  orbit  out  of  range) 


Omnidirectional  fluence/yaar  (d/em*2/MeV)t 


Ave  Energy 

Quiet 

Active 

1.9  Mev 

1 .611+0013 

1.661+0015 

2.1  MaV 

6.188*r:i2 

7.538+0012 

2.5  MaV 

3.33E+0012 

4.658+0012 

2.9  MaV 

1.79E+0012 

2.878+0012 

3.6  MaV 

8.451+001 1 

2.08+80012 

4.3  MeV 

3.308+0011 

1.061+0012 

5.7  MeV 

1.298+0011 

3.498+0011 

6.8  MeV 

6.038+0010 

1.768+0010 

8.5  MeV 

3.011+0010 

7.298+0010 

9.7  MeV 

2.008+0010 

3.858+0010 

10.7  MeV 

1 .90E+0010 

2.458+0010 

13.2  MeV 

1.068+0010 

1.218+0010 

16.9  Mev 

5.678+0009 

5.268+0009 

19.4  MaV 

3.058+0009 

2.898+0009 

26.3  MeV 

2.168+0009 

2.058+0009 

30.9  MaV 

1.478+0009 

1.528+0009 

36.3  MeV 

1.118+0009 

1.218+0009 

41.1  MeV 

5.888+0008 

6.728+0008 

47.0  MeV 

4.088+0008 

4.598+0008 

99.0  MaV 

4.708+0008 

5.128+0008 

65.7  NeV 

4.678+0008 

4.39e+0008 

81.3  MaV 

3.218+0008 

2.978+0008 

Integral  Omnidirectional  Fluence/year  <«/cm*2)  to  90.4  NeV 


Energy 

Quiet 

Active 

1.1  Mev 

1.858+0013 

2.268+0013 

1 .9  Mev 

5.688+0012 

9.32C+0012 

2.3  MeV 

3.208+0012 

6.318+0012 

2.7  MeV 

1.878+0012 

4.458+0012 

3.1  MeV 

1.198+0012 

3.30E+0012 

5.5  MeV 

3.538+0011 

7.518+0011 

5.9  Mev 

3.038+0010 

6.118+0011 

7.7  Mev 

1.958+0011 

2.94E+0011 

9.3  MeV 

1.468+0011 

1.788+0011 

10.1  MeV 

1.30E+0011 

1.47E+0011 

11.3  MeV 

1.128+0011 

1.178+0011 

15.1  Mev 

7.22E+0010 

7.13E+0010 

23.7  MeV 

4.608+0010 

4. <58+0010 

20.9  MeV 

3.47E+0010 

3.58E+0010 
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32.9  MV 

40.2  MV 

43.2  MV 
SO.B  MV 

39.2  MV 

72.2  MV 


2.892*0010 
2.078*0010 
1 .902*0010 
1.592*0010 
1.192*0010 
$.242*0009 


2.982*0010 

2.092*0010 

1.892*0010 

1.342*0010 

1.118*0010 

5.402*0009 


NASA  AP8MAX  Integral  Ormldlrtetlonal  Fluanee  for  >  90.4  MV  proton* 
7.392*0009  protona/cm*2 


32 


